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Recently, we have derived the leading and subleading isospin-breaking three-nucleon forces using 
the method of unitary transformation. In the present work we extend this analysis and consider the 
corresponding two-nucleon forces using the same approach. Certain contributions to the isospin- 
violating one- and two-pion-exchange potential have already been discussed by various groups 
within the effective field theory framework. Our findings agree with the previously obtained results. 
In addition, we present the expressions for the subleading charge-symmetry-breaking two-pion 
exchange potential which were not considered before. These corrections turn out to be numerically 
important. Together with the three-nucleon force results presented in our previous work, the results 
D ' of the present study specify completely isospin-violating nuclear forces up to the order q 5 /A 5 , where 

pL^ , q (A) denotes the soft (hard) scale. 

^ ! ■ PACS numbers: 21.45.+v,21.30.-x,25.10.+s 
> 

■ I. INTRODUCTION 

o 

O . 

iy-j i The interactions between two nucleons are one of the best studied strong interaction processes, which is due to the 
i large data base of proton-proton and neutron-proton scattering experiments. As such, many fine features of the 
strong interactions and their interplay with the electromagnetic and weak interactions can be unraveled from such 
studies, provided a sufficiently accurate theoretical tool is available to match the sometimes astonishing precision 
of the data. In this paper, we are interested in the effects of isospin violation in the two-nucleon sector. Although 
isospin is an approximate symmetry of the QCD Lagrangian, the difference in the up and down quark masses combined 
with the long-range electromagnetic force leads to appreciable deviations from the isospin limit. To quantify these 
effects, it is mandatory to have a framework that consistently incorporates these various sources. Such a scheme has 
been developed in the last decade, namely chiral nuclear effective field theory. It extends the so successful chiral 
kv^j ' perturbation theory for mesons and meson-baryon systems to processes involving 2,3,4, .. . nucleons. Chiral nuclear 
effective field theory (EFT) has already been applied to study isospin-violating two-nucleon forces (2NFs), see e.g. 
d El [J IE @ So why coming back to this topic? Recently, we have derived the leading and subleading 
isospin-breaking three-nucleon forces (3NFs) using the method of unitary transformation In the present work, 
we will apply this framework to study the corresponding isospin-breaking two-nucleon forces at the same order in 
the low-momentum expansion. Thereby, we will re-derive the many interesting results already obtained in the earlier 
studies, but also work out the contributions, which have were never been studied before. To the order we are working, 
we have thus constructed the complete set of isospin-violating few-nucleon forces. 

The material in this paper is organized as follows. In Section [n] we discuss the effective Lagrangian underlying our 
calculation and also the corresponding power counting. Section lfTTl is devoted to the derivation of the isospin-violating 
2NFs, consisting of the one- and two-pion-exchange potentials and the corresponding contact interactions. The main 
novelty is the next-to-leading order two-pion-exchange potential, which is derived here in its complete form for the 
first time and also turns out to be numerically large, very similar to the isospin-conserving case. In Section llVI we 
demonstrate explicitely the consistency between the two-nucleon forces derived here and the corresponding isospin- 
violating three-nucleon forces obtained in Ref. Section Ivl contains the summary and outlook. 
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II. POWER COUNTING AND EFFECTIVE LAGRANGIAN 



Within the Standard Model, isospin violation has its origin in the different masses of the up and down quarks and 
the electromagnetic interactions. At low energy, isospin-breaking effects in few-nucleon systems can be studied in the 
systematic and model-independent framework of chiral effective field theory. This method is based on the most general 
(approximatively) chiral invariant Lagrangian for pions and nuclcons which includes all possible interactions consistent 
with the isospin violation in the underlying theory. Consider first isospin-breaking in the strong interactions. The 
QCD quark mass term can be expressed in the two-flavor case as: 

£22? = -5?(mu + md)(l + e73)g, where e= m " ~ mrf „ _1 . ( 2 .1) 

2 m u + m d 3 

The above numerical estimation is based on the light quark mass values utilizing a modified MS subtraction scheme at 
a renormalization scale of f GeV [Uj ■ The isoscalar term in Eq. (|2.1(l breaks chiral but preserves isospin symmetry, ft 
leads to the nonvanishing pion mass, M 2 = (m u + rrid)B ^ 0, where B is a low-energy constant (LEC) that describes 
the strength of the bilinear light quark condensate. Further, this term generates a string of chiral-symmetry-breaking 
interactions in the effective Lagrangian which are proportional to positive powers of M 2 . The isovector term (oc T3) 
in Eq. 1)2. 1|1 breaks isospin symmetry and generates a series of isospin-breaking effective interactions oc (eM 2 )" with 
n > 1. It is, therefore, natural to count strong isospin violation in terms of eM 2 . 1 

Electro mag netic terms in the effective Lagrangian can be generated using the method of external sources, see e.g. 
[Til IT2T for more details. All such terms are proportional to the nucleoli charge matrix Q = e (1 + t^)/2, where 
e denotes the electric charge. 2 More precisely, the vertices which contain (do not contain) the photon fields are 
proportional to Q n (Q 2n ), where n — 1,2,.... Since we are interested here in nucleon-nucleon scattering in the 
absence of external fields, so that no photon can leave a Feynman diagram, it is convenient to introduce the small 
parameter e 2 ~ 1/10 for isospin-violating effects caused by the electromagnetic interactions. 

In the present study we adopt the same power counting rules for isospin-breaking contributions as specified in |9( . In 
particular, we count 

1 4 

£ ~ e ~A' (4^~A^' (2 - 2) 

where q (A) refers to a generic low-momentum scale (the pertinent hard scale). The iV-nucleon force receives 
contributions of the order ~ (q/k) v , where 

v = -4 + 2n 7 + 2N + 2L + ^ ViK . (2.3) 

i 

Here, L and Vi refer to the number of loops and vertices of type i and n 7 is the number of virtual photons. Further, 
the vertex dimension Aj is given by 

Aj = di + in, - 2 , (2.4) 

where m is the number of nucleon field operators and di is the g-power of the vertex, which accounts for the number 
of derivatives and insertions of pion mass, e and e/(47r) according to Eq. i|2.2f> . Finally, we adopt the counting rule 
q/m ~ (9/A) 2 for the nucleon mass m, which ensures that all iterations of the leading-order NN potential contribute 
to the scattering amplitude at leading order (q/A)° and thus have to be resumed, see for more details. 

Let us now specify the terms in the effective Lagrangian we will need in the present work. It is given in terms 
of the nucleon isodoublet and the isovector pion field tv. Utilizing the heavy baryon framework, the relevant 



1 Notice that isospin-breaking effects are in general much smaller than indicated by the numerical value of e, because the relevant scale 
for isospin— conserving contributions is the chiral-symmetry-breaking scale A x rather than m u + . 

2 Or equivalently, one can use the quark charge matrix e (1/3 + Tg)/2. 
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isospin-symmetric terms in the effective Lagrangian in the nucleon rest-frame are [Hi IT?| : 



£(0) 



ido + ^ra ■ Vtt - Jt_r ■ (tt x tt) 



A^ + 



<5m + 4 Cl M 2 -^M 2 7r 2 ' ° 2 ^ 2 ' C3 



F 2 



2F 



2 ^-ijk ^abc 

CiT a (Vj 7Tb) (Vfc 7T C ) 



iV- 



£(2) = N i 



V »gA 

2m 4mi 7 ' 



2d 



16 - »18 „ , 2 



tct • Vtt + 8id 28 M 2 d Q 



N - 



(2.5) 



where M denotes the pion mass to leading order in quark masses, M = B(m u + md), and -F can be identified with 
the pion decay constant in the chiral limit and with the electromagnetic interactions being switched off. Further, m 
denotes the physical value of the average nucleon mass, m — (m p + m n )/2, which is related to the bare mass m via 

o 

77i =77i +8m. Notice that we are using the physical and not the bare nucleon mass in the heavy baryon expansion, 
see ^| for more details. This leads to the unusual <5m-term in Eq. I|2.5|l . In addition, qa denotes the axial-vector 
coupling constant and c,; and dj (dj) are further low-energy constants (LECs). The relevant isospin-violating part of 
the Lagrangian reads 0, |2(], |2J : 



£(2) = ~^C(n 2 -7r 2 3 ) + Nl[2c 5 eM 2 T 3 -^eM 2 (Tr-T)TT 3 



N ■ 



£(3) = N i 
£(4) = N t 



he\*l - tt 2 ) + ^e 2 F 2 r 3 + e 2 ((7r • t)tt 3 - 7r 2 r 3 ) 



2d\-j — dig — 2d lg n 2 



eAra ■ Vtt 3 



N 



2ff3 + 94 



e 2 F B ■ Vtt 3 + ^-e 2 Fa- V vr 3 r 3 + ig 13 e 2 F 2 (l + r 3 )d Q + ^eM 2 a ■ V [t x tt] 3 



4e 



28 



8e 39 eM 4 r 3 



TV- 



(2.6) 



where C, /i, and ej are further LECs. Here, several comments are in order. First, we do not include in Eqs. H2.6(l 
the ess,A0~ and (714^5-terms which do not lead to isospin-breaking vertices with no pion fields. Secondly, for the sake 
of simplicity, we refrain from showing terms with four pion fields in the Lagrangians in Eqs. I|2.5|l . I|2.6I) . The explicit 
form of such terms is of no relevance for our work. We will, however, briefly discuss these terms in section Till Al In 
addition, we do not show in Eqs. I|2.5[) . (|2.ti|) isospin-violating NN contact interactions, which will be discussed in 
detail in section ITlI Dl Notice further that the complete form of the Lagrangian has not yet been worked out. 
We will discuss the relevant structures from in section 1111 Al Finally, we do not consider terms with photon 
fields. Those terms give rise to long-range electromagnetic interactions between two nucleons, which are extensively 
discussed in Refs. (23, |U I24I |2*H and will not be considered in the present work. Notice that the effects of these 
interactions are enhanced at low energy due to their long range, see e.g. |2(| for more details. In addition to these 
purely electromagnetic forces, -^-exchange contributions have to be taken into account at the order v — 4. The 
explicit expressions for the corresponding NN potential can be found in 



III. ISOSPIN-BREAKING NN FORCE 



There are many ways to derive nuclear forces from the effective Lagrangian in Eqs. (|2.5|l . I|2.6|l . In this work, we will 
use the method of unitary transformation [2jfl which leads to energy-independent and hermitean potentials. Let us 
first briefly remind the reader on the main idea of this approach. A system of an arbitrary number of interacting 
pions and nucleons can be completely described by the Schrodinger equation 

H\*)=E\V), (3.1) 

where H denotes the Hamilton operator, which specifies the interaction of pions and nucleons and can be obtained 
from the Lagrangian using the canonical formalism. Notice that due to the creation of pion field quanta via terms in 
H, the state might contain components with an arbitrary number of pions. Instead of solving the above infinite - 
dimensional equation for few-nucleon system, it is advantageous to project it onto a subspace of the Fock space, 
that contains only nucleonic states. The resulting equation can be solved using the standard methods of few-body 
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physics. Let rj and A be projection operators on the states \(j>) and which satisfy rj 2 = r/, X 2 = A, r/X = A^ = and 
A + i] = 1. Eq. (|3.1|) can then be written in the form 



rjHrj rjHX \ / 



E 2) ■ 0.2) 



We are now looking for the unitary operator U which has to be chosen in such a way that the transformed Hamilton 
operator is block-diagonal: 



In [27j we have adopted the following ansatz for the operator U 



(3.3) 



_ / „(1 + AtA)-Va _At(l + 4At)-i/a x 
^ ~ ^ A(l + AU)- x /2 A(l + A4t)-V2 ; > ^ 

which goes back to the work by Okubo 28]. The operator A in the above equation has only mixed non- vanishing 
matrix elements: A = XArj. We stress that the parameterization of the unitary operator U in Eq. I|3.4|l is not the 
most general one. The effective Hamilton operator acting on the purely nucleonic subspace of the Fock space can 
then be obtained via 

H cS = rjHrj = 77(1 + A^A)~ 1/2 (H + A+ H + HA + X s HA) (1 + A t A)- 1/2 r, , (3.5) 

The requirement in Eq. (|3.3I) leads to a set of coupled equations for the operator A, which can be solved perturbatively 
within the low-momentum expansion along the lines of |27j . One then ends up with a set of operators that contribute 
to -ffeff at a given order in the low-momentum expansion. These operators are constructed out of vertices in the 
effective Lagrangian and corresponding energy denominators, see 

[HHll for more details. The expressions for the 
nuclear potential are obtained by evaluating two-nucleon (2N), three-nucleon (3N) etc. matrix elements of these 
operators. Let us now be more specific and consider various isospin- violating contributions to the nuclear force up to 
order v = 5. 



A. One— pion— exchange potential 



The isospin-conserving one-pion-exchange (1PE) potential has been studied to one loop in |l8| using both the S- 
matrix approach, which relies on the standard technique used in quantum field theoretical calculations, and the 
method of unitary transformation. In that work, we restricted ourselves to isospin-invariant contributions. We now 
extend this analysis and include isospin-violating corrections to the 1PE potential up to v — 5. While the main focus 
of |l8fl was to study the quark mass dependence of the nuclear force, here we are only interested in the physically 
relevant case and do not need to consider the chiral expansion of the various LECs. Thus, there is no need to evaluate 
explicitly all loop diagrams which lead to pion and nucleon mass and wave-function renormalization as well as to 
renormalization of the pion-nucleon vertices. In the following, we will explain how to perform the complete calculation 
of the 1PE potential including renormalization of various LECs within the method of unitary transformation and work 
out the general structure of the isospin-breaking 1PE potential up to the considered order. 

First, we introduce, similar to |18| . renormalized pion fields and masses in the following way: 

tt* = Z; ± 1/2 tt ± , z„± = i + SZ^ , M 2 ± = M 2 + SM 2 , 

7T r ° = Zj /2 tt° , Z 7V o = l + SZ n + bZ^ , M 2 = M 2 + SM 2 - 5M 2 . (3.6) 

The quantities 5Z^ and 5M 2 denote isospin-invariant contributions to the pion wave-function and M 2 while SZ^ and 
SM 2 represent the corresponding isospin-breaking terms. These quantities can be expanded in powers of the generic 
low-momentum parameters as follows: 



SZ n = SZW + SZjp + . . . , 5Z„ = SZW + . . . , 

SM 2 = {SM 2 )^ + {SM 2 )^ + ... , SM 2 = (SM 2 )^ + {SM 2 )^ + ... , (3.7) 
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(a) (b) 



(c) 



(d) (e) (f) 



FIG. 1: Various contributions to the pion mass and wave function renormalization. Graphs (a) and (b) are 
the leading isospin-invariant contributions, while the diagrams (c) and (d-f ) show the leading and subleading 
contributions involving insertions of isospin-violating vertices. Dashed lines refer to pions; solid dots and 
filled diamonds represent isospin-invariant vertices with Aj =0 and 2 while crossed circles and crossed 
diamonds denote isospin-breaking vertices of dimension Aj —2 and 4, respectively. 



where superscripts correspond to the power of the small parameters according to Eq. (|2.2() . The leading isospin- 
invariant corrections 62^ and (dM^)^ result from graphs (a) and (b) in Fig. [I] and have already been considered 
within the method of unitary transformation in |l8j| . Leading and subleading isospin-violating contributions are given 
by diagrams (c) and (d)-(f) in Fig. ^ respectively. The leading contribution to the charged-to-neutral pion mass 
difference is entirely of electromagnetic origin and given by the C-term in Eq. (|2.6() : 

(SM 2 J^ = ^C. (3.8) 

Here and in what follows, = 92.4 MeV refers to the measured value of the (charged) pion decay constant. 3 Notice 
that graph (c) only contribute to the charged-to-neutral pion mass shift and that there are further pion self-energy 
corrections due to virtual photons, which are not shown explicitly in Fig. ^ The experimentally known pion mass 
difference M w ± - M w o = 4.6 MeV allows to fix the value of the LEC C, C = 5.9 • 1CT 5 GeV 4 . Notice that the natural 
scale for this LEC is F 2 A 2 /(47r) 2 ~ 3 • 10 -5 GeV 4 if one adopts A — M p . We do not need to explicitly evaluate 
higher-order corrections to the pion mass and wave function renormalization given by diagrams (d)-(f) in Fig.^ The 
relevant contributions are incorporated by using physical values for the charged and neutral pion masses, and the 
single-pion Hamilton operator expressed in terms of pion creation and destruction operators a\ and Oj has the usual 
form: 

H o =E/ (0r°l(*)°*(*) \I& + M h (3-9) 

where = M v ±, M3 = M„o. Notice that at the order considered the effects of the isospin-violating pion wave- 
function renormalization only shows up via an additional interaction 

H' 4 ' =-2±^Ni n ff-V* 3 N, (3.10) 

which arises from the gA -vertex being expressed in terms of renormalized pion fields and has the same structure as 
the (74 -term in Eq. Ij2.6|l . Here and in what follows, we will always work with renormalized pion fields and therefore 
omit the superscript r. 

1/2 1 /2 

Similar to the pion fields, one can define renormalized proton and neutron fields via N p = Z p N p , N„ = Z n N n . 

In the isospin symmetric case, the leading contribution to Zn results from pion loop and the d28 _ term in Eq. I|2.5|) 
(for a detailed discussion of wave function renormalization in the heavy baryon approach, see Refs. |29 | .|3C | ). Clearly, 
the NN potential derived using the method of unitary transformation includes contributions from renormalization of 



3 The difference between the charged and neutral pion decay constants is (F ± — F o)/F ~ (g/A) 4 , see e.g. Ilfl . Isospin-breaking effects 
due to F ± ^ F^o in the 1PE potential can be accounted for by small shifts in the pion— nucleon coupling constants as explained below. 
The corresponding corrections to the 2PE potential enters at order v = 6 and will not be considered in the present work. 
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external nucleon lines, which, therefore, do not need to be considered separately. Further, we remind the reader that 
the isospin-invariant nucleon mass shift 5m receives contributions at various orders in the low-momentum expansion: 

5m = 5m {1) + 5m {2) + . . . , 8m {l) ~ O (^JT^j ■ ( 3 - n ) 

The leading contribution, 5m^ = —4ciM 2 , is clearly due to the second term in the second line of Eq. while the 
subleading one, 8m^ , receives contributions from pion loops as well as from the counterterms proportional to LECs 
/i,2,3- In addition to isospin-invariant shifts, there are also isospin-breaking shifts Sm to the nucleon mass m^- 

m N ^( m * ^ ) = m + Urn r 3 . (3.12) 

The leading and subleading contributions to the proton-to-neutron mass difference are of the strong and electromag- 
netic origin, respectively: 

$ m (V = -Ac 5 eM 2 , <5m (3) = -f 2 e 2 F 2 . (3.13) 

At the order we are working, the LECs C5 and ji can be fixed from the strong and electromagnetic shifts to the 
nucleon mass: 

(m p - m n ) stI = (5m) stI = -2.05 ± 0.3 MeV , 

(rrip - m n ) cm = {5m) cm = 0.76 ± 0.3 MeV , (3.14) 



which leads to llz 



c 5 = -0.09 ±0.01 GeV- 1 , / 2 = -0.45 ± 0.19 GeV" 1 . (3.15) 



The values for the strong and electromagnetic nucleon mass shifts are taken from |3lj| . The electromagnetic shift is 
based on an evaluation of the Cottingham sum rule. In principle, this contribution could also be evaluated in chiral 
perturbation theory including virtual photons. While the formalism exists (see e.g there are still some 

subtleties to be addressed Therefore, we consider the electromagnetic mass shifts for the ground state baryon 
octet collected in |3l| the best values available. We stress that there are further corrections to 5m at higher orders 
due to pion loop diagrams and counter term insertions. We, however, do not need to consider such higher-order 
corrections in the present work since we are interested in isospin-violating corrections to the NN potential up to order 
v = 5. Indeed, since the isospin-invariant 2PE potential starts to contribute at v = 2, the corrections resulting from 
insertions of the <5m' 2 '- and 8m!^ -vertices start to contribute at orders v = 4 and v — 5, respectively. This can 
easily be verified using Eq. (|2.3(l . As will be shown below, the corrections to the 1PE potential due to the proton- 
to-neutron mass difference are either proportional to 5fhjm or to (5m) 2 . Consequently, Srh^ first contributes to 
the 1PE potential at order v = 6 which is beyond the scope of the present work. To the order we are working, the 
single-nucleon Hamilton operator takes the form: 



s,t=±l/2 J 1 ^ 



r —2 

h t dm 

2m 



(3.16) 



where n\ t (p) (n s ,t(p)) refers to the creation (destruction) of a nucleon with the spin and isospin quantum numbers 
s and t and momentum p, respectively. We stress that further isospin-breaking corrections oc 5m jm of kinematical 
origin have to be taken into account in the single-nucleon Hamilton operator entering the Schrodinger equation. For 
our purposes, however, the expression for Hq in Eq. (|3.16|) is perfectly sufficient. Notice further that here and in 
what follows, we will not separate the leading and subleading contributions 5fh^ and 5m^ to the nucleon mass shift: 
5fh ~ 5fh^ + 5m( 3 \ 

Before going into discussion of various isospin-violating contributions to the 1PE potential, it is instructive to recall its 
general (i.e. without assuming the isospin limit) form based on the phcnomenological pseudovector (PV) Lagrangian 
£pv which, for instance, in the case of the neutral pion coupled to protons has the form 



£pv = VH P ,o (N^ l5 N p ) . (3.17) 
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(a) 



(b) 



(c) 



FIG. 2: Leading (y — 2) and subleading (y = 3) isospin-violating contributions to the 1PE potential. 
Solid lines refer to nucleons and the crossed rectangle denotes an isospin-violating vertex of dimension 
Aj = 3. A light-shaded circle inserted at a pion or a nucleon line refers to a single insertion of SM% or 
5m' 2 ' + <5m < - 3 - 1 , correspondingly. Diagrams which result from the interchange of the nucleon lines and from 
the application of the time reversal operation are not shown. For remaining notation see Fig. 



Here fpp^a is the corresponding pseudovector coupling constant. Similarly, one can define the coupling constants 
fnn-n , fpmr+ and fnp-n- which correspond to the neutral pion coupled to neutrons and charged pions coupled to 
nucleons. In the case of exact isospin symmetry these are related to each other via: 

--J-f -1 

Notice that the charged pion mass enters Eq. I |3.17f t just as a scaling factor in order to make / dimensionless. The 
1PE potential can then be expressed in a general form as: 



fp 



fn 



; fnp7T 



(3.18) 



V ln (nn) = flV{M^), 

V lv {np) = -f%V(M n o) + (-l) I+1 2f*V(M„±), 



(3.19) 



where we have introduced the constants / 2 = / pp7r o / pp7r o , / 2 = / nn7r o / nn7r o , / 2 = -/p p7r o/„ n7r o and 2/ c 2 = 
— fnpir- fpn-K+ ■ Further, I = 0, 1 denotes the total isospin of the two-nucleon system and V(M i ) is defined as: 



V(Mi) 



47r ((?i • q)(a 2 ■ q) 



Mi ± q 2 + Mf 



(3.20) 



where we use the static approximation for nucleons and neglect all relativistic corrections. Charge symmetry implies 
the same interaction in the pp and nn case, i.e. fpp-K = — fnrm or fp = fn = fo ■ In case °f charge independence, one 
has fp = / 2 = /o = / 2 = f 2 . The coupling constant / 2 is related to the nucleon axial vector coupling gA via 



J 4tt 



9aM„± 



2F W 



(3.21) 



where S denotes an isospin-conserving Goldberger-Treiman discrepancy. 4 In the general case when isospin symmetry 
is not conserved, we can introduce in a close analogy to Ref. the quantities S p , S n and S c corresponding to / 2 , / 2 
and / 2 , respectively. It is sufficient to know the values of these three constants 5 P , S n and S c in order to determine 
completely the expressions for the 1PE potential in Eqs. I|3.19|l . since the coupling constant /q can be expressed as: 



fo = iz 



1 

4tt 



gAMn± 



(i + * P ) (i + <y n ) 



(3.22) 



Notice that the dominant contribution to the Goldberger-Treiman discrepancy is generated by the dis-term in 
Eq. H2.5f) and does not beak isospin symmetry: 



9A 



(3.23) 



4 The Goldberger-Treiman discrepancy is defined in terms of pseudoscalar coupling constant g. Pseudoscalar and pseudovector couplings 
lead to the same expression for the 1PE potential on-energy-shell provided the coupling constants are related via / = gM n ± /(2m). 
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In what follows, we will use the convenient form of the 1PE potential given in Eqs. (|3.19() which already incorporates 
the dominant isospin-breaking effects due to the pion-mass difference and charge dependence of the pion-nucleon 
coupling constant. 

We are now in the position to discuss the isospin-violating 1PE potential. The leading and subleading contributions 
at orders v — 2 and v — 3 are shown in Fig. [21 Since we use Sfh and are not separating 8m^ and Sfh^ , graph (b) 
in Fig. |21 contains both the order v = 2 and v = 3 contributions to the 1PE potential. Using the method of unitary 
transformation introduced above and utilizing the notation of Refs. [j| |3^|, one finds the following result for the 
diagrams (a) and (b) in Fig. 



H (0) 



A 1 



(H — E v ) 



A 1 



Hq — E,y ) 



H (0) 



(3.24) 



where r\ and r\' denote the projectors on the purely nucleonic subspace of the Fock space, while A z refers to the 
projector on the states with i pions. Further, is the leading isospin-invariant ttNN vertex proportional to gA, 
H denotes the single-particle Hamilton operator, H = Hq +H N , and E v (E„i) refers to the energy of the nucleons in 
the state 77 (77'). The contribution of graph (a) to the 1PE potential can be obtained by evaluating the corresponding 
matrix element of the operator in Eq. 1|3.24|) in the limits m — > 00 and Sfh — > 0. Clearly, this contribution is 
already included in Eqs. (|3.19|l . Similarly, the contribution of the diagram (b) is obtained by evaluating the matrix 
element in the limits m — ► 00 and 8M% — ► and expanding energy denominators in powers of Sfh. The term linear in 
Sfh leads to vanishing matrix elements, so that there is no contribution from graph (b) to the 1PE potential. Finally, 
the contribution of the last diagram (c) can be obtained by evaluating the matrix elements of the operator 



-2 V 



H (0). 



A 1 



(Hq — E v 



A 1 



(H — ) 



T] + h. c. , 



(3.25) 



in the limits m — > 00, 5m — > and bMl -> 0. Here, H^ denotes an isospin-violating vertex from the Lagrangian 
£(3) j n (|2.(j|) which is proportional to the combination of the LECs Idyj — d\% — 2di 9 . It provides a contribution 
to the quantities 5 P and 5 n : 



= -5( 3 ) = 2 2dl7 ~ rfl8 



2d 



9A 



19 eM 2 



(3.26) 



which leads to the charge symmetry breaking 1PE potential. 



Isospin-violating contributions to the 1PE potential at order v = 4 are depicted in Fig. [3J Notice that pion tadpole 
graphs with the Weinberg-Tomozawa vertex lead to vanishing contributions and are not shown in Fig. In addition, 
pion loop diagrams (but not of the tadpole type) with one insertion of the Weinberg-Tomozawa or cs-vertex do not 
contribute since only odd functions of the loop momentum enter the corresponding integrals. We do not show those 
diagrams in Fig. [3] either. Finally, we also ignore one-loop nucleon self-energy contributions with insertions of the 
pion mass difference, since they do not lead to isospin-violating contributions. Let us begin with the first graph (a) 
in Fig- 01 The contribution proportional to di$ leads to renormalization of the nucleon axial vector coupling constant 
gA while the one proportional to dis provides a dominant contribution to the isospin-conserving Goldberger-Treiman 
discrepancy, see Eq. Ij3.26|l . Further, the leading relativistic I/m-correction vanishes as does the corresponding 
isospin-invariant contribution. This can easily be understood in terms of Feynman diagrams. Indeed, the 1/m-vertex 
in the last line of Eq. (|2.5|l contains a time derivative of the pion field. Since the four-momentum is conserved, it 
gives a contribution proportional to the nucleon kinetic energy which is of the order v = 6. For graph (b) we find that 
the contributions proportional to d\§ and ciis vanish. The leading relativistic correction oc Sfh/m corresponding to 
this diagram can be obtained by evaluating matrix elements of the operators in Eq. I|3.25|) . where corresponds 
now to the second term in the last line of Eq. (|2.5|l . Taking the limit m — * 00 in Eq. I|3.16|l . expanding the energy 
denominators in powers of Sfh and keeping only terms linear in Sfh, we find: 



v ™ - 1 2^ \Sk 



9A 



Ti X T2J3 



1 



(<j\ ■ qi)(a 2 ■ (p2 +P2')) + (<?i ' (Pi +Pi')){&2 ■ qi) 



(3.27) 



where f>i (pi ') denotes the incoming (outgoing) momentum of the nucleon i and <fi = p\ — p± = —(pi — P2). The 
order Aj = 2 isospin-invariant interaction in the next two graphs (c) and (d) is due to the d28 _ vertex, Sm^ 2 ' and 
the nucleon kinetic energy. The <i28 - term in the Lagrangian is proportional to the nucleon equation of motion and is 
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00 



(c) 



(d) 



(e) 



(f) 



t 
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00 



(i) 



G) 



00 



(1) 



4 v 



(m) 



00 



(o) 



(P) 



0) 



00 



FIG. 3: Order i/ = 4 contributions to the 1PE potential which contain isospin-violating vertices. For 
diagrams with insertions of the nucleon mass shift only one representative topology is depicted. Further 
graphs with nucleon mass insertions at different internal or external nucleon lines are not shown. For 
remaining notation, see Figs. 00 



only needed for renormalization purpose. 5 In the method of unitary transformation, this term is eliminated from the 
Lagrangian by an appropriate field redefinition, see [l8| for more details. The contributions from graphs (c) and (d) 
proportional to the nucleon kinetic energy can be obtained evaluating matrix elements of the operators in Eq. i|3.24[l , 
expanding in 6M%, 1/m and Sfh and keeping only terms proportional to SM 2 /m and Sfh/m, respectively. We find 
that graph (c) leads to a vanishing result while graph (d) provides the following contribution to the 1PE potential: 



^r=-*2^2prj [rixra], { ,, + M , y (Pi "J* "Pi +^ ), (3-28) 

Introducing the total momentum P of the two-nucleon system, P = pi + P2 = pi +p~2 , the above expression can be 
re-written as: 



^ = ' — ( ttf- \ [t> ^ r,-; '- 1 : • />. (8.2!)) 



Thus, this potential vanishes in the two-nucleon center-of-mass (CMS). The contributions from graphs (c) and 
(d) proportional to 8mS 2 > are found to vanish. As pointed out before, the contribution from diagram (e) in Fig. 
proportional to (SM 2 ) 2 is already taken into account in Eqs. (|oM9|) . For graph (f) we obtain the following contribution 
to the 1PE potential: 



*l - -IX-.* I MA I ,3, 31 Pi-9)p2-?) 



5 It absorbs the ultraviolet divergence in the nucleon Z-factor calculated to one loop. 



10 



Notice that the isospin-violating piece has the same structure as the correction due to the pion mass difference at 
order v = 2 but is 5M%/(5fh) 2 ~ 660 times weaker. Next, graph (g) is expected to provide a correction proportional to 
6fh6M%. We find that the corresponding contribution vanishes. Further, the contribution of diagram (h) is included 
in the 1PE potential in Eqs. I|3.19[) . where one has to account for shifts in the pion-nucleon coupling constants: 

Sf^ = e^Fl + ^ , */« = - * e^Fl + ^ . (3.31) 

9A 2 g A 2 

Consider now diagram (i) in Fig. |3| The isospin-violating counterterms of dimension Aj = 4 include both the strong 
term cx e39 which leads to the nucleon mass shift, and the electromagnetic one oc g%3 which is proportional to the 
nucleon equation of motion and absorbs the ultraviolet divergence in the nucleon ^-factor. We find that this diagram 
as well as all remaining diagrams (j)-(s) in Fig. y| either lead to vanishing contributions or renormalizc various LECs. 
In particular, contributions of graphs (k) and (i), (j), (1) can be expressed in terms of isospin-conserving and isospin- 
violating nucleon self-energy corrections, while diagrams (m), (p), (r) and (n), (o), (s) give rise to isospin-invariant 
and isospin-breaking renormalization of the pion-nucleon coupling constants. None of the pion loop diagrams lead to 
any form-factor-like behavior, i.e. have a non-trivial dependence on the momentum transfer between two nucleons. 
Thus the contribution of all these diagrams is taken into account by using the general expressions for the one-pion 
exchange potential in Eqs. I|3.19ll expressed in terms of renormalized quantities. Stated differently, the contributions 
of these diagrams only lead to charge-dependent shifts in the strength of the 1PE potential. We stress that since the 
corresponding LECs g\ and g\ are not known experimentally and have to be determined from the data, and because 
we are not interested in the quark-mass dependence of the strength of the 1PE potential, we do not need to evaluate 
the loop diagrams in Fig. [3] explicitly. Finally, we have also verified the finding of Ref. Q, that the contributions of 
diagrams (1) and (n) cancel. 

Let us now discuss the corrections to the 1PE potential at order v = 5 represented by the diagrams depicted in Fig. 01 
Notice that we again refrain from showing various kinds of diagrams which lead to vanishing contributions as explained 
above. For the same reason, we also do not show graphs with an insertion of the ititNN vertices proportional to 
/i; fi (but not of the tadpole type) as well as pion tadpole diagrams proportional to Cj with an insertion of the 
5M 2 -vertex at the pion line, which do not lead to isospin-violating contributions to the potential. First, we note 
that the contribution of graph (a) proportional to the LECs d\$ and dxs is already included in Eqs. (|3.19|) . Further, 
similar to the case of diagram (a) in Fig.yJ the corresponding 1/m-correction vanishes at this order. Diagrams (b) 
and (c) also do not lead to new structures in the 1PE potential: the term oc e?28 provides a shift of the strength 
of the charge-symmetry-breaking (CSB) 1PE potential via the nucleon wave-function renormalization while the 
contributions proportional to the nucleon kinetic energy and to 5m^ vanish. Similarly, diagrams of the type (c) and 
(d) in Fig. |3| but proportional to Sm^ instead of 5m^ 2 \ which are not shown explicitly in Fig.QJ lead to a vanishing 
result. The contribution of diagram (d) is already included in the expression for the 1PE potential in Eqs. I|3.19|l . 
Graph (e) leads to a vanishing result. Consider now the contribution of diagram (f). Order A, = 5 isospin-violating 
counter terms have not yet been worked out. Here we simply list all possible structures consistent with the usual 
symmetry constraints: 

£(«) = Aft ^ a e 2 . V [ T x n] 3 + b x eM* a ■ V V 2 ^ + b 2 e M%a ■ W 3 + ^3 eM* a ■ V7r 3 ) N , (3.32) 

where a and bi are LECs. The first term in the above expression does not contribute at order v — 5 due to the 
presence of the time derivative. The 1PE potential proportional to 6, can be cast into the form of Eqs. I|3.19|l plus 
additional CSB short-range interactions. Further, one has to take into account relativistic 1/m-terms 

£(5) = . 2d 17 - rf 18 -2d 19 2 ^ ^ _ . ^ 
zmfn 

which, however, contribute at higher orders due to the presence of the time derivative. Graph (g) in Fig. ^represents 
the correction to the 1PE potential due to an insertion of the order Aj = 5 isospin-violating counter terms. To 
the best of our knowledge, the latter have not yet been worked out. They may include eM^-terms proportional to 
the nucleon equation of motion which contribute to the nucleon Z-factor as well as terms oc M 2 e 2 /(4-7r) 2 which give 
further corrections to the nucleon mass difference. Again, we do not need to evaluate explicitly the contributions of 
this diagram as well as all remaining graphs (h)-(v) in Fig.0]since they do not lead to any new structures in the 1PE 
potential. We have verified that diagrams of type (k) only lead to shifts in the pion-nucleon coupling constants Si. 
We further note that the loop integral in graph (j) has no logarithmic ultraviolet divergence, which is consistent with 
the fact that diagram (f) has no contribution due to counter terms of electromagnetic origin. 

To summarize, isospin-violating corrections to the 1PE potential up to v = 5 are accounted for by using the expression 
in Eqs. (j3~ll?|) and further corrections in Eqs. (|3~T7|) . and PD) . The CSB corrections in Eqs. lj!H7|l and (|3~25|) 
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FIG. 4: Order v = 5 contributions to the 1PE potential which contain isospin-violating vertices. Filled 
rectangles refer to isospin-invariant vertices of dimension Aj = 1 while triangles denote isospin-breaking 
interactions of dimension Ai = 5. For remaining notation, see Figs. Q |2 E3 



obtained within the method of unitary transformation agree with the ones found in |8( using a completely different 
framework. The correction cx (Sfh) 2 in Eq. I|3.30[l can be found e.g. in [34j|, see also discussion in 



B. Two— pion— exchange potential 



Let us now discuss the leading and subleading isospin-violating two-pion-exchange potential which can be expressed 
in momentum space as: 

= (t x 3 r|) [V c + V s (8 1 ■ 8%) + V T {8x ■ q)(8 2 ■ q)] + (rf + r|) [w c + W s (8i ■ 8 2 ) + W T {dx ■ q)(8 2 ■ q)] , (3.34) 

with the six functions Vc(q), ■ ■ ■ , Wx{q) depending on the momentum transfer q = \q\. The subscripts refer to the 
central (C), spin-spin (S) and tensor (T) components in the potential. Further, Vi and Wi correspond to charge- 
symmetry conserving and charge-symmetry breaking pieces, respectively. The dominant contributions arise at order 
v = 4 from diagrams shown in Fig. |SJ In this figure, graphs (a)-(d) represent the effects due to the pion mass 
difference, graphs (e)-(h) provide contributions proportional to the nucleon mass difference, and the last two graphs 
(i) and (j) are due to a single insertion of the isospin-breaking nirNN-vertex proportional to the LEC C5. 
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FIG. 5: Leading isospin-breaking correction to the 2PE potential at order v = 4. For notation, 
see Figs. [T] Hand 01 



The isospin-violating two-pion-exchange (2PE) potential due to the pion mass difference has been considered in [|| . 
As shown in this reference, it can be expressed in terms of the corresponding isospin-invariant contributions without 
performing any additional calculations. To that aim, one can first decompose the isospin-invariant 2PE potential into 
the isoscalar and isovector pieces 

v 2n = v 2 \ + v 2 \ n • t 2 . (3.35) 

The leading isospin-breaking effects die to M n ± ^ M n o are incorporated properly if one uses M n , defined as 

= ^M w± + iil^o , (3.36) 

in the scalar part V 2 % and expresses the vector part as following: 

{V 2 „ (M„± ) for pp and nn , 
V 2 l(M„o) for np, T = l , (3.37) 
Vl(M n ) for np, T = 0. 

These results are valid modulo (<5Af 2 /M 2 ) 2 -corrections. Notice that terms oc (5M 2 /M 2 ) 2 start to contribute to the 
2PE potential at order v = 6 and thus need not be considered in the present work. Equivalently, we can write the 
2PE potential in the form 

V 2 „ - V 2n (M n ) + V 2n (M w ) Tl . T2 + — — Tl T . 2 + Ol J , (3.38) 



where M n is the average pion mass 



M v = \{M n ± +M v o) . (3.39) 



Substituting the expressions for the isospin-invariant 2PE potential at v = 2 given e.g. in [3^| into Eq. I|3.38|l we 
obtain for the non-polynomial parts of the contribution of diagrams (a)-(d) in Fig. [SJ 

(4) = SM 2 1 f 4 2 _ ( 4M 2(g ff 4 _ 4 2 2 (n 4 _ g 2 _ _ J^lMl 

c 128tt 2 F4 4M 2 + q 2 \ yA * \ tV yA yA ' y v yA yA ' 4M 2 + q 2 
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where the loop function L A (q) reads 



rA , , CO . K 2 lu 2 +q 2 s 2 + 2kqujs 

L (q = — In „, . „ ^ 

W 2q AM 2 (A 2 + q 2 ) 



uj = ^q 2 + AM 2 , 



VA 2 - 4M 2 ■ 



(3.41) 



The above expression for L A (q) is given in the spectral-function rcgularization (SFR) framework with A being the 
corresponding cut-off. The limit A — * oo corresponds to dimensional regularization (DR). 

Consider now diagrams (e)-(h) in Fig. [S] which include one insertion of the nucleon mass shift. The contributions 
of graphs (e) and (f) can be obtained within the method of unitary transformation by evaluating the corresponding 
matrix elements of the operator 



5e,5f 



2 IT 



X 1 



2 (Ho — E,y 



■H^fjH^ 



A 1 



8 (Ho — E v > 



-H<Vf}H<® 



(Hq — Efj)(H — E v > 
A 1 



. H (o) 



8 (Ho — E v >)(Ho — Efj) 



1 



X 1 X 2 
-H®* -,H^ ] -^ , 

2 (Ho — E v ) (Ho — E v ) 



(Ho — Efj)(H — E n ) 
A 1 



H (0) 



H (0). 



(Ho — E v ) 



H (o) 



T) + h. c. 



(3.42) 



in the limits m — > oo, 5M 2 — > by expanding the denominators in powers of 5m and keeping only terms linear in dm. 6 
Similarly, contributions of diagrams (g) and (h) can be obtained by evaluating the matrix elements of the operators 



V, 



5g _ 



2 V 



H (o) 



A 1 



- H (0)_ 



A 2 



(Ho - E v ) (Ho — E n ) 
A 1 A 1 



H (0) 



. H (0). 



(Ho - E v ) (H — E v ) 



h (°) + hW 

H (o) 



- H (0)_ 



X 1 



(Ho — E n ) (Hq — E n 
T) + h. c. , 



. H (0) 



and 



V 5h - -W 



ff (0). 



A 2 



(Hq — E n i 



. H (0) + ff (0) 



A 2 



(Ho 



_ H (o) 



V, 



respectively. Finally, the operators corresponding to the last two graphs in Fig. read 
Vg=rf 



H (o) ?l H (o) a2 h (2) , H m 



LU (UJI + LU 2 ) 



X 2 X 1 X 1 X 1 



(ll>i + LU 2 ) LU 
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to 



and 



v 2 5 i = v' 



A 2 



(U>1 + LU 2 ) 



_ H (o) +H (o)_ 



X 2 



(iOi + UJ 2 )' 



(3.43) 



(3.44) 



(3.45) 



(3.46) 



where the w's refer to the pionic free energy in the isospin limit and H^ 2 ' denotes the isospin-violating ttttNN vertex 
from the Lagrangian £^ proportional to the LEC C5. Performing a straightforward evaluation of the matrix elements 
of these operators, we obtain the following result for the leading CSB 2PE potential: 



W, 



(4) 



9 a \l9\5mMl 



it: 



(4) 



GAttF* { AM 2 + q 2 

9 A $ m 
327rFi 



Ag 2 A 6rh-(5m)^ (2M 2 + q 2 )A A (q) 



1 w (4) - 9 * Sffl A A (a) 



(3.47) 



Notice that the contribution oc SM^ discussed above can also be obtained from Eq. 13.421 if one takes the limits m — > 00, 8m — > in 
H . 
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FIG. 6: Subleading isospin-breaking correction to the 2PE potential at order v = 5. For notation 
see Figs. EHl 



where the loop function A A (q) is given by 

A A tq) = — arctan — . (3.48) 

Here, several comments are in order. First, one should keep in mind that we again only show explicitly the non- 
polynomial terms. Secondly, we found that the planar box graph (e) and the "football" diagrams (h) and (j) lead 
to vanishing contributions. Notice further that we have also included the order v = 5 contribution from graph (f) 
in Fig. The latter is proportional to the LEC f% or, cquivalently, to (<5m) em and has the same structure as the 
contribution of diagram (i) in Fig. [5] which is oc C5 or, equivalently, oc (Sfh) stI . It is therefore convenient to combine 
these contributions and to express them in terms of e.g. 5fh and (<5m) str . Finally, we would like to emphasize that the 
result given in Eq. (|3.47|l agrees with previous calculations, see 0, EH and for related older work. 

Consider now the subleading isospin-breaking two-pion-exchange potential generated by diagrams shown in Fig. [S] 
First, charge-symmetry conserving contributions from graphs (a) and (b) due to the pion mass difference can be 
obtained using Eq. (|3.38[) from the corresponding isospin-invariant 2PE potential. We find: 

V^=~VP=-^ a \c,A\q). (3.49) 



The contributions of graphs (c) and (d) result in the method of unitary transformation from the operators 



v 2tt 



V 

2 ' 



H (o) 



A 1 



A 1 



(Ho — E v ) (H — E v ) 
A 1 „/i-s A 1 



tf(O) 



(Ho — E v ) (Hq — E v ) 



. H (0) 



(Ho — E. q ) (Hq — E rj 
7} + h. c. , 



. H (0) 



and 



vl d = -J 



H [ 



\ 2 \ 2 



(Hq — E r y ) 



(H — E n ) 



. H (o) 



f] + h. c. 



(3.50) 



(3.51) 



where H^> refers to ttttNN vertices from £W proportional to Cj. Further, contributions from graphs (e) and (g) 
in Fig. can be obtained from V 2 * in Eq. (|3.46|) by replacing H^ by H^ and H ^ by H ^ , respectively. The 
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contribution of diagram (f) results from V 2 5 ^ in Eq. I|3.45[l with being replaced by and has already been 
included in Eqs. (|3.47|) . Finally, contributions of graphs (h), (i) and diagram (j) can be obtained by evaluating matrix 
elements of the operators V 2 °^' 5 in Eq. (|3.42|) and in Eq. (|3.43|) . respectively, where one of the vertices is 
replaced by the isospin-breaking ttNN vertex from Eq. I|2.6|) and all possible permutations of the operators 
and are taken into account. Notice that in this case one only needs to keep the pionic free energy in 

the corresponding denominators. We found that diagrams (d), (g) and (j) lead to vanishing contributions while the 
subleading CSB potential generated by diagrams (c), (e), (h) and (i) in Fig. reads: 



W, 



(5) 



1 



967r 2 Fi 



L\q){ -g\5fh 



48A^(2ci+c 3 ) 
4M2 + o 2 



w4 5) = 



-Ws ] 
q 2 b 



9 2 a 



167T 2 Fi 



g\ 8m (18ci + 2c 2 — 3c 3 ) 
18c 3 ) - (28m 



g\ 8m (5c 2 
L A (q) Urn a 



(28m - (5m) str ) (6ci - c 2 - 3c 3 ) 
(Sm) str ) (ca + 6c 3 ) 



5A, 



(3.52) 



where /3 = eM 2 (2di7 — ^18 — 2dig). Notice that we have also included the contribution of diagram (e) in Fig. 
with the C5~vertex being replaced by the /2-vertex from Eqs. I|2.6[) . which appears formally at order v = 6. 7 This 
contribution has precisely the same structure as the one proportional to the LEC C5 with the overall strength 2(8fh) cra 
instead of (6m) str . In Eqs. (|3.52ll we have expressed 2(Sm) cnl + (Sih) stI as 28m — (Sm) str . 

To summarize, the charge symmetry conserving 2PE potential is due to the pion mass difference and includes a 
central component at order v = 4 and tensor and spin-spin components at order v = 5 given in Eqs. I|3.40|l and 
(|3.49|l . respectively. The CSB 2PE potential has all central, tensor and spin-spin components at orders v = 4 and 
v = 5, and the corresponding expressions are given in Eqs. I|3.47|l and l|3.52|l . Further, we stress that our results for 
the isospin-violating 2PE potential are consistent with taking 



An 



2F\ 



(3.53) 



in the isospin-conserving 2PE potential. This expression already accounts for the Goldberger-Treiman discrepancy. 
Finally, to the order we are working, the constant /3 can be expressed in terms of the pion-nucleon coupling constants 
as follows: 



f 2 - f 2 

J p J n 



(3.54) 



C. Two— pion— exchange potential in coordinate space 



Let us now take a look at the isospin-violating potential in coordinate space. Similar to Eq. I|3.34|l we define: 



V 2 rr = (r x 3 r 2 3 ) V c + V s (ai • ff 2 ) + V T S 12 + (7? + t£ ) W c + W s {&i ■ S 2 ) + W T S 12 



(3.55) 



where Si 2 = 3 a\ • f a 2 ■ r — ai • a 2 and the functions Vc(r), ■ ■ ■ , Wt(t) depend on the distance r. Consider now the 
charge-symmetry-conserving 2PE potentials V$ (r) which are due to the pion mass difference and correspond to , 



and Vg defined in Eqs. I|3.4U|I and l|3.49|l . In order to obtain the r-space expressions (at r ^ 0), one usually 
switches to the spectral-function representation for the non-polynomial part of the potential which has the form: 

Im[Vi(-ifi)] 



2 r 
VM = - / 



dpi [i ■ 



/i 2 + q A 



(3.56) 



The corresponding contribution oc f\ does not lead to isospin breaking and is not considered. 
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Here, Im[Vi(— i/i)] is the mass spectrum (or the spectral function) entering this representation, which results from 
the analytical continuation of the momentum-space functions to q = + — Using this representation one obtains 
for the functions Vi(r) for r > 0: 



V c (r) = 



27T 2 7 



2M„ 



dmie-^lm\y c {-i^)] , 



Vr(r) = - 



V 8 (r) = - 



1 



g 7r 2 r 3 
1 



2 A / 



dfi n e~^ r (3 + 3^r + ^r 2 ) lm[V T (-in}] 



6ir r 



f dune-"* (n 2 Im[Vr(-«At)] - 3 Im^H/i)] ) . 

■/ 2M_ V 7 



(3.57) 
(3.58) 
(3.59) 



For the CSC 2PE potential oc 5M 2 the spectral-function representation takes a slightly more complicated form as 
compared to Eq. I|3.56|l . namely: 



Vi(q) = lim 



Im [Vi (-*/*)] 8M X Im [V, (-2*M W - ie)] 



2M»+e 



M 2 + <? 2 



4M 2 + q 2 



(3.60) 



The second term in this expression accounts for the corresponding shifts of the threshold when pions of different 
masses are exchanged between the nucleons. Notice that both terms in Eq. (|3.6U|I go to infinity for e — > 0, but their 
sum is, of course, finite in this limit. Although one can, in principle, use this representation on order to obtain Vi(r), 
we choose different strategy by using Eq. (|3.38l) in coordinate space. Substituting the expressions for the leading and 
subleading isospin-invariant 2PE potential in r-space j^S in Eq. (|3.38|l , we find: 



V#>(r) 



A^oo 



vf\r) 



5M 2 M„ 



256tt 3 F4 r 2 
5M 2 g 2 A c 4 e- 2 * 



[4<&(-l + g\)xK (2x) + (-1 - 6g 2 A + g\(U + ix^K^x)] , 



192tt 2 .F 4 r 4 
5M 2 g 2 A c A e~ 2 * 



" '■l-' + -'-'-"-^^5 : (8 + y(5 + y)) 



967T 2 ^i 



-(1 + 2s(l + x)) + ^§4? ^ 2 + ^ 2 + »)) = 
r 4 1927r z i< 4 r 4 



(3.61) 



where x = M„-r, y = Kr and If* are the modified Bessel functions. In the case of , we could only perform the 
integral in Eq. (|3.57|) analytically for A — > oo which corresponds to the DR result. 

To obtain the CSB 2PE at orders v = 4 and v = 5 in r-space, it is convenient to use Eqs. (|3.57() - (|3.59(l (with Vi being 
replaced by Wi). In that case, the spectral-function representation in Eq. (|3.56|) is valid. Using Eqs. 13.47|) . (|3.52|) 
we find: 



g\ e 2x 



({5m) str (l + xf - 25m g\{2 + x{A + x{2 + x)))) 



iuf(r) 



9 a 



5127r 2 Fi r 4 



{{5m) st * - A5mg 2 A ) (2 - 2x 2 + y{2 + y)) , 



5m g\ e 2x t , ,v 5m g\ < 



384tt 2 F 4 r 4 
5m g A e~ 2x 



(4 + x(5 + 2x)) + 



(l + 2x + 2x 2 )) 



768tt 2 F 4 r 
Sfhg A e~ y 



_(8 + y(5 + y)) 



384tt 2 F 4 r 4 



(2 + 2y + y 2 )) 



(3.62) 



and 



A^oo 



g\ 5mx(-5c 2 + 18c 3 + 4(2ci + c 3 )x 2 ) 



32tt 3 F 4 r 4 
+ (25m - (<5m) str ) x (c 2 + 6c 3 )J^o(2.t) 

+ (g\ 5m (-5c 2 + 18c 3 + 2(6c x - c 2 + 5c 3 )x 2 



17 




_jgO I i i i i I i i i i I i i i i I i i i i I i i i i I _g I i i i i I i i i i I i i i i I i i i i 

1 1.2 1.4 1.6 1.8 2 2 2.5 3 3.5 4 

r [fm] r [fm] 



FIG. 7: Central (upper row), tensor (middle row) and spin-spin (bottom row) components of the 
charge-symmetry-conserving (CSC) 2PE potential in coordinate space, utilizing DR (dashed) and 
SFR (solid lines). The SFR result corresponds to A = 700 MeV. Left /right panels: Distances from 
1...2/2...4 fm. 



Wf\r) 



A^oo 



A^oo 



+ (25m - (5m) stT ) (c 2 + 6c 3 + 2(2 Cl + c 3 )x 2 )) Kx(2x) 



9a m * 
'96TT 3 F£r 4 

9a m * 



(c 4 8m + (3 g A ) (l2x K {2x) + (15 + Ax 2 )^ [2x] 



{c i 5rh + /3g A )[SxK (2x) + (3 + 2x^)K 1 (2x) 



(3.63) 



Again, at order !/ = 5we could only perform integrals in Eqs. (|3.57|l - (|3.59|l analytically for A — > oo. 

It is now interesting to compare the strength of the corresponding r-space potentials. Here and in what follows, 
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FIG. 8: Central (upper row), tensor (middle row) and spin-spin (bottom row) components of the 
charge-symmetry-breaking 2PE potential in coordinate space. Long-/short-dashed lines: DR at 
v = 4 and v = 5; Solid/dot-dashed lines: SFR at v = 4 and v — 5. The SFR results correspond to 
A = 700 MeV. 



we adopt the same values for the LECs Cj as in our work [2(j: ci = —0.81 GeV 1 , C2 = 3.28 GeV 1 , C3 = —3.40 
GeV" 1 and c 4 = -3.40 GeV" 1 . Further, g A = 1. 27, M OT = 138.03 MeV and = 92.4 MeV. For the strong nucleon 
mass shift (<5m) str , we use the value given in Eq. 1)3.14(1 . Finally, in our numerical estimations we set (3 = since the 
value of this LEC is not known at present. Notice, however, that a 1% relative deviation between f p and /„ leads 
to ~ 6ffiC4, so that the strength of the resulting CSB potential is comparable to the one of the CSB potential 
oc C4. The CSC and CSB 2PE potentials are plotted in Figs. [7| and |SJ respectively, for two choices of the cut-off A in 
the SFR: A = 700 MeV and A = 00 which is equivalent to DR. The results for V^\r) and W^\r) for A = 700 MeV 
are obtained numerically. First of all, we notice that all CSC 2PE contributions have similar strength ~ 100 keV at 
r = 1 fm and ~ 4 — 6 keV at r = 2 fm. Although the central potential Vq^{t) is formally dominant, the subleading 

contributions V^ s {r) are numerically large due to the large value of the LEC C4. To get further insights into the 
importance of the CSC 2PE potential, it is instructive to compare its strength with the strength of the corresponding 
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FIG. 9: Ratio of the CSC tensor 27r-exchange (using SFR with A = 700 MeV) and l7r-exchange 
potentials as a function of r. 



1PE potential resulting due to the pion mass difference, which provides the dominant isospin-violating contribution 
to the NN force at order v = 2. It has the following form in momentum space: 

^ (<?) = - (, 27V ) 6M " + Mfp ■ (3 ' 64) 
The corresponding r-space expressions read (for r > 0) : 



e 



In Fig. |51 we have plotted the ratio /Vffl as a function of r. The 2PE contribution is significant for r < 2 fm. 
At larger distances it becomes negligible (less than 1% for r > 3.5 fm) compared to the 1PE contribution due to its 
shorter range. 

Let us now switch to the CSB contributions displayed in Fig. |S] Similar to the CSC potential, the subleading terms 
oc Cj are numerically enhanced due to large values of these LECs. The strongest contribution is given by the subleading 
central potential which reaches ^ 150 — 300 keV at r = 1 fm and ~ 6 keV at r = 2 fm and is dominated by terms 
oc C3. This is similar to the isospin-symmetric case, where the subleading central 2PE potential is known to be very 
strong. To enable a more detailed comparison between the leading and subleading contributions to the CSB 2PE 
potential, we plot in Fig. llOl the corresponding ratios. While the subleading central component is significantly stronger 
than the leading one in a wide range of distances r, the tensor and spin-spin components are of the same size as the 
leading ones only for r < 2 fm. One should, however, keep in mind that our results obtained within the heavy baryon 
formalism become formally invalid at very large distances. This problem with the heavy baryon formalism has been 
first observed in the single-nucleon sector and can be dealt with using e.g. the Lorentz invariant scheme proposed 
by Becher and Leutwyler [3jJ, see also |4jj. It is clear, however, that the NN interaction due to two-pion exchange 
becomes very weak at large distances, so that the problem with the formal inconsistency of the heavy baryon approach 
is expected to have little relevance for practical applications, see [^J for more details. 

Finally, we would like to emphasize that the SFR results for CSC and the leading CSB contributions are rather close 
to the ones obtained in DR even for 1 fm < r < 1.5 fm, see Fig. [7] Stated differently, the corresponding potentials 
seem to be strongly dominated by long-range components in 27r-exchange. Deviations between the SFR and DR 
results at short distances are significantly larger for CSB contributions at order v = 5. This is consistent with the fact 
that new CSB short-range terms with two derivatives start to contribute at this order. One expects that such new 
counterterms will largely remove the dependence of observables on the cut-off in the spectral-function representation. 
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FIG. 10: Ratio of the subleading (y = 5) and leading (y — 4) CSB 2PE potential as a function 
of r. Solid/dashed/dot-dashed line: central/tensor/spin-spin component. All results correspond to 
SFR with A = 700 MeV. 



D. Contact terms 



We now discuss short-range isospin-breaking interactions. The leading contact terms in the NN potential at v 
are proportional to the quark mass difference and violate charge symmetry: 



K ( o2t = 03,! (rf + rf ) + fa,* <?i • S 2 (rf + r 2 3 ) , 



(3.66) 



where 03,-j ~ £ M 2 /(F 2 A 2 ) are constants. We remind the reader that the subscripts of the spin and isospin matrices 
refer to nucleon labels while the superscripts denote the corresponding vector indices. Both terms in Eq. I|3.66[l lead to 
the same structure in the potential when anti-symmetrization with respect to the nucleons is performed. Consequently, 
it is sufficient to keep only one term. 8 At order v = 4 one has to take into account isospin-violating short-range 
interactions without derivatives of electromagnetic origin. In addition to the terms which have the same structure as 



the ones in Eq. 



and thus only provide order ~ e 2 / (Air F^) 2 shifts of 03, i, new CSC interactions appear: 



04,2 <?i • a 2 rf r| 



(3.67) 



Again, one of these two terms can be eliminated performing anti-symmetrization of the potential. Finally, at order 
v = 5 one needs to include CSB terms with two derivatives which are proportional to the quark mass difference. In 
the two nucleon CMS these terms read: 



Vfit = (rf+rf) 



05,1 q 2 + 05,2 k 2 + (05,3 q 2 + 05,4 k 2 )(ai • a 2 ) + -05,5 (<?i + a 2 ) ■ (k X q) 



+ 05,6 (<?i • q)(a 2 ■ q) +05, 7 (^i ' k)(a 2 ■ k) 



+ «05,8 (rf - r|) k X q ■ {&x - a 2 ) 



(3.68) 



Equivalcntly, one can apply a Fierz transformation in the corresponding Lagrangian in order to eliminate redundant terms. 
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+ i0 5 , 9 [ti x -r 2 ] 3 k x q - [a x x ct 2 ] + 05,io [ti x t 2 ] 3 ({a x ■ q)(a 2 ■ k ) - (a x ■ k)(a 2 ■ qj) , 

where 0$^ ~ eM^/(F%k A ) are further constants and fc = (p+p')/2. Performing anti-symmetrization of this potential, 
it is easy to see that half of the terms proportional to 05, x , /3s, 2, 05, 3, 05,4, 05,6 and 05, 7 are redundant. Similarly, 
terms proportional to 05,&, 05,9 and 05,w, which lead to mixing between the T = 1 and T = states, generate the 
same structure when the potential is anti-symmetrized. We are, therefore, left with five independent short-range 
terms at order v — 5 (for example, one can take terms proportional to 05, 1, 05,3, 05,5, 05,6 and 05, s). 

The isospin-breaking short-range terms up to order v = 5 feed into the matrix-elements of the S- and P-waves in 
the following way: 

C So, PP \Kontl 1 So, PP) = P 1 P SO+01SO(P 2 +P 2 ) , 

CSo, nn ItUtlX nn) = 0[" s l o ~ 0iso (p 2 + p' 2 ) , 

( 3 P , pp\V cont \ 3 P , pp) = -( 3 P , nn\V cont \ 3 P , nn) = /3 3PO pp' , 

( 3 Pl, pp|Kont| 3 Pi, PP) = ~( 3 Pi, nn\V cont \ 3 P x , nn)=0 3P1 pp' , 

( 3 Pi, Pp\V cont \ 3 P 2 , PP) = -( 3 P 2 , nn\V cont \ 3 P 2 , nn)=0 3P2 pp' , 

( 1 P 1 , np\V cont \ 3 P x , np) = /3ipi-3pi PP' , (3.69) 

where the new LECs /3^ , /3™g , /3isoj 03PO, 03P1 and /3ipi_3pi can be expressed in terms of linear combinations 
of the LECs 03,i, 04 t i and 05^. Notice that we have adopted here the convention according to which the np matrix 
elements (with exception of the last term in Eq. (|3.t)9l) ) do not change by switching off isospin violating contact terms. 

As pointed out in 5], one should, in principle, also take into account isospin-violating CSC contact terms associated 
with the contributions to the 2PE potential cx M%. These terms arise since the derivative in Eq. (|3.38|) has to be 
applied not only to the non-polynomial part of V 2 \{M^), but also to the corresponding counterterms which depend on 
M^. The resulting CSC contact terms have fixed coefficients (in terms of qa and F^) and are of the order ~ 6M%/M% 
instead of expected ~ SM^/A 2 compared to the corresponding isospin-invariant terms. We have found that these 
contact interactions cancel against the ones which result from taking the derivative of the non-polynomial part of 
V^^Mrr). Thus, Eqs. I|3.40[l and l|3.49[) give the complete result for the corresponding CSC 2PE potential, and no 
additional contact terms need to be taken into account. Finally, we stress that no contact interactions depending on 
the total two-nucleon momentum P (similar to Eqs. I|3.27fl . (|3.28l) ) appear up to the order v = 5. 



IV. CONSISTENCY OF THE 2N AND 3N FORCES 



It is well known that two- and three-nucleon forces have to be consistent with each other. In [£| and |43 we have 
derived isospin-violating 3NF at orders v = 4 and v = 5. It has 1PE, 2PE and contact pieces and results from 
insertions of pion and nucleon mass differences as well as the C5- and /i j2 -vertices in Eq. (|2.tj|) . Two 3NF diagrams 
out of seven shown in Fig. 1 of Ref. |42| appear to be of a special interest due to the fact that they include reducible 
topologies, i.e. time-ordered graphs with purely nucleonic intermediate states. These two diagrams are depicted in 
Fig. Illland lead to the following contributions to the 3NF at order v = 4: 



+ Qi ■ Qj [(Ti ■ T k )r 3 - (-n • Tj)tI] I 



= J2^™Ct(§t) ( g/+^ 2 fa x ^ ft x ft] • * ■ (4-70) 



where i, j and k are nucleon labels and q\ = Pi' — Pi- Further, CV is one of the two leading order four-nucleon LECs 
[l5| . Notice that these nonvanishing contributions are in strong contrast to the corresponding isospin-invariant 3NF 
forces oc oi and oc g\ Ci as well as isospin-breaking ones oc 6M% g\ and oc 8M 2 g\ C{ which are known to vanish, 
see e.g. ^J, 0, 0. It is, therefore, instructive to look at the origin of these particular 3NFs in more detail. To 
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(a) 



(b) 



FIG. 11: Isospin violating 3NF diagrams which include reducible topologies. For notation, see 

Figs, mil 



that aim we evaluate the 3N scattering amplitude oc Sm g\ and cx Sm g\ Ct and compare the result with the iterated 
2N potential. We will demonstrate that the iterated 2NF + 3NF given in Eq. (|4.7U|) reproduces correctly the 3N 
scattering amplitude. This is an excellent and rather non-trivial check of our results. 

The most convenient way to evaluate the scattering amplitude is using the Feynman graph technique. In the following, 
we will use the method suggested in Ref. 8], in which the proton-to-neutron mass difference is removed from the 
Lagrangian by an appropriate field redefinition in favor of new isospin-violating terms oc Sm which are easier to 
handle in practical applications. Only two of such new isospin-breaking terms may lead to contributions oc Smg\ 
and oc Smg 2 A Ct in the amplitude: 

c — 

C = Sm(7r x n) 3 - — NUa-p, (t x n) 3 }N . (4.71) 
4iv m 

Notice that in contrast to nuclear forces, the on-shell scattering amplitude we are interested in is unique and does 
not depend on field redefinitions. The relevant Feynman diagrams are shown in Fig. 1121 Notice that Feynman graphs 
resulting from diagrams (a) and (b) with one of the g^-vertices being replaced by the vertex corresponding to the 
second term in Eq. I|2.5|l lead to contributions oc 1/m which are irrelevant for our discussion. The contribution from 
graph (a) is given by: 



T 



12a 



gA \ i 2 i 2mA_| 



. 2F„ J [q, 2 - M 2 + ie] 2 [p k 2 - m 2 + ie] [q 3 2 - M% + ie] 

x (<?i ' qi)(-Vk ■ qi){&k ■ qj){-°j ■ qj) (-25m) {qi) t abi t\tI (r k ■ Tj) , (4.72) 

where we use the relativistic nucleon propagator. Here, A + is the projection matrix onto positive-energy states and 

Pk —p'k + Qi = P'k+Pi~Pi = C\J Pk' 2 + m2 + y ^ Pi' 2 + m2 ~ V Pi 2 + m2 i Pk' + q-^j is the four-momentum of the nucleon 
k in the intermediate state. Notice further that the amplitude is multiplied by i in order to match with the standard 
normalization of the nonrelativistic T-matrix. The consistency of this procedure can be verified e.g. by calculating 
the static isospin-invariant 1PE potential. Using the relation: 

i2mA+ ' -oflV (4.73) 



\pk 2 — m 2 + ie] 5T + ie \m 
where ST = (pi' 2 + Pk' 2 — pi 2 — Pk 2 ) /2m, we can express T 12a as (modulo l/m 2 -corrections): 

Pi Pi \ 1 



T 12a = ^2Smv 2 k [T. l XT k ] s (^-^r) vl^Tk-Tj), (4.74) 



where 



2m 2m / ST + ie 



9 A \ 2 (gj • qi ){d k ■ qj ) 
2Fj {q 2 + M 2 Y 



(4.75) 



We now use the equality 



(pi' 2 


A 2 ) 




\ 2m 


2m j 





P? P k ' 2 , Pu 2 , 5T \ (4 ?6) 



2m 2m 2m 2m 
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FIG. 12: Feynman diagrams contributing to the CSB 3N scattering amplitude. Graphs resulting 
from the interchange of the nucleon lines and/or application of the time reversal operation are not 
shown. Empty circles denote insertions of the 5m-vertices from Eq. (14. 711 . 



to rewrite Eq. (|4.74|) in the form: 



T 



12a 



5m 



2 -T 2 - /2 

-Pk -Pi 



i 5m vf k [n x r fc ] 3 vlj (x fc • Tj) 



■Pk 



5T 



~ V kj ( T k ' T 3 



(4.77) 



The first term in the above equation can be identified with the iteration of the isospin-violating 1PE potential in 
Eq. (|3.28|) between the nucleons i and k and the isospin-invariant static 1PE potential between the nucleons k and 
i) Vin — vh (T k • tj). The second term in Eq. I|4.77[) thus corresponds to the genuine contribution of the 3NF. 
Performing the algebra for spin and isospin matrices, we obtain the following result: 



i 5mvf k [n x T k f vL (t* • Tj) = 5m 



{di ■ g^jdj ■ cfj) 



9A 

2FJ (^ + M2)2(^.2 + M2) 

Wi X ^f] ■ a k [Ti X Tj} 3 + (fr ■ Qj 



(4.78) 



(Ti ■ T k )rf - ( Ti ■ Tj)4 



j] 3 + i[q l x qj] ■ a k (r, • T k )r 3 - (t, ■ Tj)r^ 



To get the complete expression for the 3NF we need to take into account the contribution of the diagram resulting 
from graph (a) in Fig. 1121 by interchanging the ordering of the pion propagators and summing over the nucleon labels. 
This leads to cancellation of the terms in the third line of Eq. I|4.78|l . The final result agrees with V lla in Eq. 14.70fl . 

Similarly, the scattering amplitude corresponding to the graph (b) in Fig. 1121 is: 



T 



lib 



-i c L5mv\ k [n x T k 



I 3 — 
1 \ 2m 



Pi 



1 



2m 5T 



C's + C T {oj ■ <7 fc ) 



(4.79) 



Notice that we use the following Feynman rule for contact terms: (— i) [Cs + Ct (<Jj • <?&)], which leads to the 2N 
potential V — Cs + Ct (5j ■ a k )- We rewrite the amplitude T 12b in the form: 



T 



12/. 



. 5m o 3 
1 2m~ V * k [Tl X Tk] 



2 ~ 2 
Pi - Pk 



-.,2 ,2 
Pi +Pk 



5T 



C s + C T [Zj ■ dvO 



i5mv1 k [Ti x T k } 3 



C's + C T (Bj ■ ffk) 



(4.80) 



Again, the term in the first line of the above equation gives the iterative contribution to the amplitude, while the 
term in the second line is the genuine 3NF contribution: 



i5mvf k [Ti x T k ] 3 



Cs + C T ((Jj ■ S k ) 



5m 



9A 

2F„ 



[Ti X Tfc 



(<fi 2 + M2)2 

ji C s (ak -qi)+i C T (a 3 ■ $) - C T [aj x a k ] ■ gij 



(4.81) 
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The first two terms in the curly brackets cancel against the contribution of the time-reversed diagram, so that the 
final result agrees with V llb in Eq. I|4.70|l . 

Finally, it is easy to see that the amplitude corresponding to graphs (c) and (d) in Fig. ^] is reproduced by the 
iteration of the isospin-breaking 1PE potential in Eq. Ij3.27|l and the leading isospin-symmetric 2N potential with no 
need for an additional 3NF. We thus have verified the consistency of our results for isospin-violating 1PE 2N force 
and 3NFs. 



V. SUMMARY AND OUTLOOK 



The pertinent findings of this study can be summarized as follows: 

1) We have applied the method of unitary transformation to study the isospin-breaking NN forces. We have 
derived all forces (in the absence of virtual photons) up to v = 5 in the chiral expansion. Together with the 
3NFs worked out in Ref. 0, this completely specifies the isospin-breaking few-nucleon forces up to this order. 

2) The isospin-violating one-pion exchange potential is a combination of terms of order v — 2,3,4 and 5. It is 
specified in Eqs. I|3.19|l . I|3.27|l . (|3.28|l and (I3.3U|) . In particular, we have reproduced the expressions for the 1PE 
potential, class IV, found recently in Ref. jg. 

3) The results for the isospin-violating 2PE potential in momentum space are given in Eqs. H3.40JI . i|3.47Jl . (|3.49(1 
and 13.5211 . The corresponding expressions in coordinate space are discussed in detail in section Rll CJI In 
particular, we have reproduced the expressions for the 2PE potential at order v — 4 worked out in Rcfs. 7, 36]. 
The main new result is the 2PE potential at order v = 5 which has not been considered in EFT before. We 
find that the derived subleading CSB 2PE at order v = 5 gives numerically the dominant contribution to the 
CSB TPE potential, cf. Figs. 181101 As in the isospin-conserving case, this effect can be traced back to the large 
magnitude of the dimension two pion-nucleon LECs 02,03 and 04 (which is well understood, see |46f1. 

4) The contact interactions appear at orders v = 3,4,5 and are listed in Eq. I|3.69[) . These are parameterized in 
terms of the LECs pi and pi. 

5) We have also demonstrated explicitely that our results for the isospin-violating 2N and 3N forces are consistent 
with each other. 

These results pave the way for new precision studies. To make this point more transparent, we reiterate that the 
finite-range part of the potential, i.e. the 1PE and 2PE pieces, depend on the mass shifts 6M% and Sin which are 
well known, and the strong nucleon mass shift (<5m) str which is less well known. For the latter quantity, one can use 
the value based on the Cottingham sum rule, see Eq. I|3.14|l . In addition, the pion-nucleon coupling constants f£, 
/„ and are not well known at present. In principle, these constants can be extracted from an independent partial 
wave analysis of the pp and np data, such as the new Nijmegen PWA 47]. Provided such an extraction is possible, 
the finite-range part of the isospin-violating nucleon force is completely determined. One can then try to fix the 
values of the LECs Pi and Pi from the low partial waves in the np and pp systems in order to make predictions for 
the corresponding nn partial waves. It remains to be seen whether this ambitious program can be carried out. Also, 
one should keep in mind that the isospin-conserving NN forces have so far only been worked out up to order v = 4 
in the power counting. 



Acknowledgments 



E.E. would like to thank Prof. W. Glockle for helpful discussions and warm hospitality during his stay in Bochum, 
where part of this work has been done. This work has been supported by the U.S. Department of Energy Contract 
No. DE-AC05-84ER40150 under which the Southeastern Universities Research Association (SURA) operates the 
Thomas Jefferson National Accelerator Facility and by the Deutsche Forschungsgemeinschaft through funds provided 



25 



to the SFB/TR 16 "Subnuclear Structure of Matter". This research is part of the EU Integrated Infrastructure 
Initiative Hadron Physics Project under contract number PJI3-CT-2004-506078. 



[1] U. L. van Kolck, Soft physics: Applications of effective chiral lagrangians to nuclear physics and quark models, PhD thesis, 

University of Texas, Austin, USA, 1993, UMI-94-01021. 
[2] U. van Kolck, J. L. Friar and T. Goldman, Phys. Lett. B371, 169 (1996), [nucl-th/9601009]. 

[3] U. van Kolck, M. C. M. Rentmeester, J. L. Friar, T. Goldman and J. J. de Swart, Phys. Rev. Lett. 80, 4386 (1998), 
[nucl-th/9710067]. 

[4] E. Epelbaum and U.-G. Meifiner, Phys. Lett. B461, 287 (1999), [nucl-th/9902042] . 

[5] M. Walzl, U.-G. Meifiner and E. Epelbaum, Nucl. Phys. A693, 663 (2001), [nucl-th/0010019]. 

[6] J. L. Friar and U. van Kolck, Phys. Rev. C60, 034006 (1999), [nucl-th/9906048]. 

[7] J. L. Friar, U. van Kolck, G. L. Payne and S. A. Coon, Phys. Rev. C68, 024003 (2003), [nucl-th/0303058] . 
[8] J. L. Friar, U. van Kolck, M. C. M. Rentmeester and R. G. E. Timmermans, Phys. Rev. C70, 044001 (2004), [nucl- 
th/0406026]. 

[9] E. Epelbaum, U.-G. Meifiner and J. E. Palomar, Phys. Rev. C71, 024001 (2005), [nucl-th/0407037]. 
[10] H. Leutwyler, Phys. Lett. B374, 163 (1996), [hep-ph/9601234]. 
[11] R. Urech, Nucl. Phys. B433, 234 (1995). 

[12] U.-G. Meifiner and S. Steininger, Phys. Lett. B419, 403 (1998), [hep-ph/9709453]. 
[13] G. Miiller and U.-G. Meifiner, Nucl. Phys. B556, 265 (1999), [hep-ph/9903375]. 
[14] S. Weinberg, Phys. Lett. B251, 288 (1990). 
[15] S. Weinberg, Nucl. Phys. B363, 3 (1991). 

[16] V. Bernard, N. Kaiser and U.-G. Meifiner, Int. J. Mod. Phys. E4, 193 (1995), [hep-ph/9501384]. 

[17] N. Fettes and U.-G. Meifiner, Nucl. Phys. A693, 693 (2001), [hep-ph/0101030]. 

[18] E. Epelbaum, U.-G. Meifiner and W. Glockle, Nucl. Phys. A714, 535 (2003), [nucl-th/0207089] . 

[19] S. Steininger, Reelle und virtuelle photonen in chiraler Storungstheorie, PhD thesis, University Bonn, Germany, 1999. 

[20] N. Fettes, U.-G. Meifiner, M. Mojzis and S. Steininger, Ann. Phys. 283, 273 (2000), [hep-ph/0001308]. 

[21] J. Gasser, M. A. Ivanov, E. Lipartia, M. Mojzis and A. Rusetsky, Eur. Phys. J. C26, 13 (2002), [hep-ph/0206068]. 

[22] E. A. Ueling, Phys. Rev. 48, 55 (1935). 

[23] L. Durand III, Phys. Rev. 108, 1597 (1957). 

[24] G. J. M. Austin and J. J. de Swart, Phys. Rev. Lett. 50, 2039 (1983). 
[25] V. G. Stoks and J. J. de Swart, Phys. Rev. C42, 1235 (1990). 

[26] E. Epelbaum, W. Glockle and U.-G. Meifiner, Nucl. Phys. A747, 362 (2005), [nucl-th/0405048]. 
[27] E. Epelbaum, W. Glockle and U.-G. Meifiner, Nucl. Phys. A637, 107 (1998), [nucl-th/9801064]. 
[28] S. Okubo, Prog. Theor. Phys. 12, 603 (1954). 

[29] G. Ecker and M. Mojzis, Phys. Lett. B410, 266 (1997), [hep-ph/9705216]. 

[30] S. Steininger, U.-G. Meifiner and N. Fettes, JHEP 09, 008 (1998), [hep-ph/9808280]. 

[31] J. Gasser and H. Leutwyler, Phys. Rept. 87, 77 (1982). 

[32] J. Gasser, A. Rusetsky and I. Scimemi, Eur. Phys. J. C32, 97 (2003), [hep-ph/0305260]. 
[33] E. Epelbaum, W. Glockle and U.-G. Meifiner, Phys. Lett. B439, 1 (1998), [nucl-th/9804005] . 

[34] V. G. Stoks, The magnetic moment interaction in NN phase shift analysis, PhD thesis, University Nijmegen, The 
Netherlands, 1990. 

[35] E. Epelbaum, W. Glockle and U.-G. Meifiner, Eur. Phys. J. A19, 125 (2004), [nucl-th/0304037]. 
[36] J. A. Niskanen, Phys. Rev. C65, 037001 (2002), [nucl-th/0108015]. 
[37] S. A. Coon and J. A. Niskanen, Phys. Rev. C53, 1154 (1996). 

[38] N. Kaiser, R. Brockmann and W. Weise, Nucl. Phys. A625, 758 (1997), [nucl-th/9706045] . 
[39] T. Becher and H. Leutwyler, Eur. Phys. J. C9, 643 (1999), [hep-ph/9901384]. 

[40] T. Fuchs, J. Gegelia, G. Japaridze and S. Scherer, Phys. Rev. D68, 056005 (2003), [hep-ph/0302117]. 

[41] R. Higa, nucl-th/0411046. 

[42] E. Epelbaum, nucl-th/0412003. 

[43] S. N. Yang and W. Glockle, Phys. Rev. C33, 1774 (1986). 

[44] J. A. Eden and M. F. Gari, Phys. Rev. C53, 1510 (1996), [nucl-th/9601025]. 

[45] U. van Kolck, Phys. Rev. C49, 2932 (1994). 

[46] V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. A615, 483 (1997), [hep-ph/9611253]. 
[47] M. C. M. Rentmeester, R. G. E. Timmermans and J. J. de Swart, nucl-th/0410042. 



